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\Q . Abstract. Accreting X-ray pulsars are among the best observed objects of X-ray astronomy 
with a rich data set of observational phenomena in the spectral and timing domain. While the 
general picture for these sources is well established, the detailed physics behind the observed 

^\ ■ phenomena are often subject of debate. 

T—i ■ We present recent observational, theoretical and modeling results for the structure and dy- 

t-H ' namics of the accretion column in these systems. Our results indicate the presence of different 
CNj ■ accretion regimes and possible explanations for observed variations of spectral features with 
' luminosity. 
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1. Introduction 

Accreting X-ray pulsars exist in binary systems, where a companion star transfers 
matter to the neutron star via Roche lobe overflow, a strong stellar wind or a Be outflow 
disk, e.g., Frank, King & Raine (2002). The sizes of these systems vary widely, generally 
in the range 10 6 to 10 s km. 

Independently of the way matter is transferred, close to the compact object the flow 
of matter is dominated by the strong magnetic field B ~ 10 12 G. The accreting plasma 
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couples to the magnetic field lines at the Alven radius, which is a few 1000 km for typical 
neutron star parameters (see, e.g., Caballero & Wilms (2012) and references therein). 
The inf ailing matter follows the magnetic field lines, forming accretion columns on the 
neutron star magnetic poles with a surface area of ^lkm 2 and a height that can range 
from meters to many km. 

The observed X-ray spectra are usually relatively simple and broadly similar across 
sources. They are normally well described by a power law with exponential cutoff, affected 
by absorption and including iron fluorescence lines which are produced within the binary 
system, as well close as far away from the X-ray source. In some accreting X-ray pulsars, 
broad absorption-like features are observed which are caused by resonant scattering of 
photons on electrons whose energies are quantized into so-called Landau levels, spaced 

by 

, eB „ B 



11.6 keV 



m c c 



10 12 G' 



For an overview of known cyclotron line sources see Caballero & Wilms (2012) 



For bright accreting X-ray pulsars Becker & Wolff (2007) have developed a model 
of the continuum spectra, where the total observed spectrum is a result of the bulk 
and thermal Comptonization of bremstrahlung, black body and cyclotron seed pho- 
tons. This model has been adapted to xspec and used for analysis of 4U 0115+63 by 



Ferrigno et al. (2009) 



Cyclotron line features are still usually fitted phenomenologically with independent 
Gaussian or Lorentzian line shapes. A more physical approach to model the complex har- 



monic cyclotron line shapes simultaneously has been pursued by Schonherr et al. (2007) 



Extensive Monte-Carlo simulations of the cyclotron resonant scattering on a large pa- 
rameter grid (following Araya & Harding 1999, 2000) have been used to create a Green's 
function based xspec convolution model, applicable to any continuum and successfully 



applied to spectra of different sources, e.g., by Suchy et al. (2008) 



In contrast to the spectra, the pulse profiles - light curves folded by the pulse pe- 
riod - can be quite unique for a given source, with the typical shape serving almost 
as a 'fingerprint'. The profiles are usually quite dependent on the energy band, with 
a common tendency to have simpler profiles at higher energies. In some sources these 
energy-dependent profiles are very stable across significant brightness variations, in other 
cases a clear evolution is visible and occasionally rather chaotic changes have been ob- 
served (e.g., Camero-Arranz et al. 2007). What drives these different modes of behaviour 
is still not well understood. 



2. Recent Results 



2.1. Pulse profile modeling 



Caballero et al. (2011) and Sasaki et al. (2012) have applied methods developed by Kraus et al. (1995) 
in order to decompose the pulse profiles of the sources A 0535+26, 4U 0115+63 and 
V 0332+53. In this approach axisymmetrical emission and the same beam pattern is as- 
sumed for both poles, but the poles can be offset from antipodal. The observed profiles 
can be explained within these models including emission from the accretion column a 
halo around the footpoint of the column as well as scattering of photons in the upper 
accretion stream. The change of emission geometry during a giant outburst of V 0332+53 



could be demonstrated by Sasaki et al. (2012) 
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2.2. Cyclotron line energy variations 
For some of the sources with cyclotron lines the derived line energies have been found to 



vary as a function of luminosity, see, e.g., Mowlavi et al. (2006) Tsygankov et al. (2007) 
Staubert et al. (2007) or Klochkov et al. (2011) The data point towards two types of 



variability: for certain sources, (e.g., V 0332+53) the centroid energy correlates negatively 
with the luminosity, for others (e.g., Her X-l) a positive correlation is observed, while 
for some sources (e.g., A 0535+26) the line appears to be at constant energies (within 
uncertainties) despite significant luminosity changes. A possible explanation has been 



developed by Becker et al. (2012) who describe different accretion regimes based on the 
definition of a critical luminosity L CIlt = 1.5 x 10 37 (£?/10 12 G) above which radiation 
pressure decelerates the matter to rest. The brightest sources ("supercritical", i.e., Lx > 
^crit), will have an accretion column where the characteristic emission height increases 
with luminosity and the observed B field decreases accordingly. In subcritical sources, 
somewhat below L cr it, Coulomb braking will be dominant and increased luminosity leads 
to a decreasing emission height and increasing observed B field. At even lower luminosities 
no significant evolution is expected. This model agrees well with the observed source 



behaviour and parameters, for details see Becker et al. (2012) 



2.3. Accretion mound stability and variations with pulse phase 

A common feature of cyclotron line sources is that the observed line energy varies with 
pulse phase. While a variation with angle is expected if one includes relativistic correc- 
tions, the observed variations are often larger than those predicted straightforwardly. 



Mukherjee & Bhattacharya (2012) and colleagues have studied the stability of accre- 



tion mounds and the distribution of the magnetic field in the mound. For mounds of a 
total accreted mass of ~ 1O~ 12 M there is an appreciable distortion of the magnetic field 
with a stronger relative field towards the rim of the mound. This distortion would have 
a strong signature in the observed spectra if mound surface were visible directly. In the 
2D MHD simulations ballooning instabilities develop which limit the maximum height 
of the mound and the amount of matter that can be contained. Beyond ~ 10 _13 M Q the 
mounds become unstable. 



3. New developments 

In order to progress significantly beyond the efforts described above, a challenging 
overall problem needs to be tackled: For each pole the X-ray continuum production 
needs to be modeled as function of energy and position in the column. This emission 
is further significantly modified by cyclotron resonant scattering strongly dependent on 
energy and angle of photon propagation. To describe the observed spectra or pulse profiles 
the emission from both poles must be combined using realistic emission geometries and 
including relativistic light bending around the neutron star. 

As one step in this process, members of our collaboration are further developing the 
Green's function models of cyclotron line scattering processes (Schonherr et al. in prep., 
Schwann et al. in prep.). Compared to previous approaches (e.g., Araya & Harding 1999, 
2000) much more general geometries are being implemented, allowing the adaption to 
light-bending models for the effectively visible emission. This work is cross-checked with 
other approaches, see, e.g., Nishimura (2008, 2011). 



In parallel an effort is ongoing to combine the modeling efforts of Becker & Wolff (2007) 



and Becker et al. (2012) with the advanced cyclotron line scattering description. In this 
'hybrid' approach, column properties and seed photons are derived from the continuum 
models, which get modified in an outer 'sheath' by the advanced cyclotron scattering 
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models. An additional component in this effort is the inclusion of ray tracing within a 
general relativistic geometry in order to describe the corresponding modifications to the 
observed spectra. 

The structure of accretion mounds is further studied in 3D MHD simulations by 
Mukherje, Bhattacharya & Mignone (in prep.). First results indicate an easy excita- 
tion of the fluting mode instability and consequent transport of matter across the field. 
This would reduce the maximum mass supported within the mounds and field distortions 
would be on small scales. 

4. Summary 

In recent years significant progress has been made on various fronts in the understand- 
ing of the physics of accretion columns in X-ray pulsars. These developments have been 
possible due to improved methods and advances in computing power. The results are 
encouraging, but also demonstrate the complexity of the issues which have to be taken 
into account. 

Within our collaboration work is now ongoing to improve the existing models and to 
combine different approaches in order to provide detailed predictions for direct compar- 
ison with the rich observational dataset. 
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